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Abstract 

The diagrammatic approach is applied to study quasiparticle transport proper- 
ties in two-dimensional d-wave superconductors with dilute nonmagnetic impurities 
both in Born and in unitary limits. It is found that a novel quantum interference 
process gives rise to a weak-antilocalization correction to the spin conductivity, 
indicating the existence of extended low-energy quasiparticle states. With comimg 
close to unitarity and the nesting, this correction is suppressed and eventually 
vanishes due to the global particle-hole symmetry. 



PACS numbers: 74.25.Fy, 73.20.Fz, 74.20.-z 

In recent years there has been increasing interest in the understanding of low- 
energy quasiparticle (QP) states in disordered d-wave superconductors The d x i_ y i- 
wave pairing state is characterized by an anisotropic energy gap, which vanishes along 
four nodal directions. There exist low-lying Dirac-type QP excitations near the gap 
nodes. In sharp contrast to the conventional s-wave superconductors, even nonmagnetic 
impurities can drastically change the behavior of the low-energy QP states in d-wave 
superconductors. A central issue, not noly experimentally relevant but also theoretically 
intricate, is whether these QP states are localized and how the disorder affects the low- 
energy QP transport properties. Over the last decade, a variety of conceptually and 
methodologically different approaches to the problem have been developed, many of 
these theories contradict each other. Based on the self-consistent treatments |^], a 
nonlinear sigma model M, or numerical studies pL pi q], some groups suggested that all 
the QP states are localized. On the other hand, Balatsky and Salkola have shown that a 
single strong impurity produces a virtual-bound state at zero energy, and the long-range 
overlaps between these impurity states yield an extended QP band [0. The singularity 
in the density of states (DOS) at zero energy obtained recently by the nonperturbative 
T-matrix method signals the QP derealization as well||. The possible appearances 
of critical states || 10 and localization-delocalization transitions []Tl| in random Dirac 



fermions have been also discussed. As a result, the problem of QP (de)localization 
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in disordered <i-wave superconductors still remains controversial and deserves further 
scrutiny. 

The problem of low-energy QP localization was usually studied or discussed on the 
basis of calculations of the DOS[||, |(| ||, ||. An alternative approach to this issue is 
to calculate the spin conductivity, which was first proposed and done by Senthil et 
al. based on the nonlinear sigma modelQ, and recently used by Zhu, Sheng and Ting 
based on the numerical transfer matrix method Q. It is well known that the quantum 
interference effects (QIEs), resulting from the cooperon and diffuson in diagrammatic 
language, play an important part in the low-energy QP transport in two-dimensional 
(2D) disordered (i-wave superconductors^, ||, 12, 113]. Unlike in a normal metal, every 



cooperon and diffuson mode in the retarded-advanced (RA) channel entails a corre- 
sponding mode in the retarded-retarded (RR) or advanced-advanced (AA) channel due 
to the local particle- hole symmetry (LPHS) in the superconducting state 0]. In the 
unitary limit and nesting case, each of these 0-mode cooperons and diffusons has a 



7r-mode counterpart induced by the global particle-hole symmetry (GPHS)|l3|, |l5fl . Re- 
cently, the above symmetries in these Goldstone modes have been used to account for 
the physical origin of previous contradictory theoretical predictions for the DOS in 
disordered d-w&ve superconductors [p^| , It is highly expected that the features of these 
Goldstone modes could be used to address the issue whether the low-energy QP states 
are localized. 

In this Letter we present a diagrammatic study of the QIEs on low-energy QP 
transport properties in 2D d-w&ve superconductors with dilute nonmagnetic impurities 
both in Born and in unitary limits. At the one-loop level, we find a new impurity- 
scattering polarization diagram related to the LPHS, which has never been considered 
in the previous diagrammatic analyses. This novel quantum interference process is 
found to have a profound effect on the spin conductivity. It is shown that, in general, 
the spin conductivity is subject to a weak-antilocalization correction while the electrical 
conductivity has a weak-localization correction. In the singlet superconductors, the 
spin of QPs is a good quantum number but the charge is not|J. Therefore, this weak- 
antilocalization effect indicates the existence of extended low-energy QP states. With 
coming close to unitarity and the nesting, the corrections of both spin and electrical 
conductivities are suppressed, and eventually vanish due to the GPHS. A semiclassical 
picture involving interfering trajectories of the novel quantum interference process is 
also presented. 

Let us start from a most extensively studied model for a 2D d x 2_ y 2-w&ve super- 
conductor, in which the normal-state dispersion and energy gap are given, respec- 
tively, by £k = —t(cosk x a + cos k y a) — fi and Ak = Ao(cosfca;a — cosk y a), with 
t the nearest-neighbour hopping integral, a the lattice constant, and \x the chemi- 
cal potential. In the vicinity of the four gap nodes k n = (=tfcp, ±/sf)/\/2, the QP 
spectrum = + A^) 1 / 2 can ^ e linearized as ek ~ [(vy-k) 2 + (v 9 -k) 2 ] 1//2 , where 
v/ = (<9£k/<9k)k n , v 9 = (<9Ak/9k)k n , and k is the momentum measured from the 
node k n . Consider pointlike nonmagnetic impurities to be randomly distributed with 
low concentration m and the impurity potential V, then the time-reversal and spin- 
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rotational symmetries are preserved (symmetry class CI[14]). In the self-consistent 
T-matrix approximation, the QP self-energy can be expressed in the Nambu spinor 
representation asfl5| ^ R ^ A \e) = riiT R ^ A \e) = (Ae =F 27)10 + r/jr^ for |e| <C 7. Here 
A is the mass renormalization factor, 7 is the impurity-induced relaxation rate, rj is 
a dimensionless parameter, and tq and Tj (i = 1,2,3) stand for the 2x2 unity and 
Pauli matrices, respectively. A use of Dyson's equation yields the impurity-averaged 
one-particle Green's functions as 



-,R(A) f ^ _ [(1 - A)e ± 27)7-0 + A k n + g k r 3 

[ (1 _ A)e±Z7 ]2_ e 2 



The impurity-induced DOS at zero energy is calculated as po = — (l/7r)Im ^ k TrG k (0) = 
4l~f/ir 2 VfV g , where I = ln(r/7) > 1 with T ~ ^/VfV g /a. The parameters 7, A, and rj 
can be evaluated consistently |l5|l , yielding 7 = 2nj/7r/7o(l + rj 2 ), A = (1 — rj 2 )(l — 
l)/(r] 2 + 21 — 1), and r\ = 2/irpolI with C7 the effective impurity potential given by 
U^ 1 = V^ 1 + X^k£k(e k + 7 2 ) _ • The Born and unitary limits correspond to rj 2 S> 2/ 
and 77 — ► 0, respectively. 

In the generic situtations, only the 0-mode cooperon and diffuson contribute to 
the QIE@, |l|. Owing to the LPHS, T 2 G R (e)T 2 = -G A (-e), they exist both in 
RA and RR channels, and can be expressed as D(q;e,e') RR ^ = C*(q; e, e') RR ^ = 
J2t C(q; e, e') RR( - A Ki n, where 

C{^e' ) r {A) =cf R{A) /[D q 2 - l {e±e% (2) 

with c RA = c RA = -c RA = c RA = -c RR = c RR = c RR = c RR = i^/irpo and 
D = (v 2 + v 2 )/&l^ the QP diffusion coefficient. Equation (2) is valid both in Born and 
in unitary limits. As in the study of disordered interacting electron systems [16], all the 



leading polarization diagrams responsible for the QIE can be generated from the lowest- 
order self-energy corrections Figs. 1(c) and 1(d) in Ref . Jl5[ . Figure 1 represents all the 
one- loop diagrams with 0-mode cooperon. Diagrams 1(a) and 1(b) have been studied 
by Altland and Zirnbauer in the context of the random matrix theory of mesoscopic 



normal/superconducting systems [ 14 1. Figure 1(c) is a new polarization diagram, its 
contribution can be shown to be of the same order as that of Figs. 1(a) and 1(b). 
It can be shown that Fig. 1(d), as well as all the one- loop diagrams with 0-mode 



diffuson, has a vanishing contribution [17]. Therefore, to the lowest order, it is Figs. 
l(a)-l(c) that result in the QIE. The nonvanishing contributions of both diagrams 1(a) 
and 1(c) stem from the existence of the cooperon in RR channel. We wish to emphasize 
here that diagram 1(c) describles a novel quantum interference process in 2D <i-wave 
superconductors. As will be shown below, it is the existence of diagram 1(c) that leads 
to a weak-antilocalization correction to the "classical" value of the spin conductivity. 

The Kubo formula is used to calculate the QP transport coefficients. For the 
spin conductivity, each vertex of the diagrams in Fig. 1 contains a vector A = 
( v g T i + v / r 3)/2[^l]- The contributions of these diagrams to the spin conductivity 
can be expressed by <r£ = (l/27r)Re(II^ yl — n RR ) (x = a ib,c), where Tl RA and IL RR 
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stand for the corresponding zero-frequency spin current-current correlation functions 



in RA and RR channels, respectively. The expressions for 14^; are given by 
n « A = \ E E C ^)u R ^ (A k G^7iG V.G^A k Gfl , 



kq 



kq i 



q kk' »j 



with G^> = G£ w (0) and C(q)^ l/lj = G(q; 0, 0)f R(A) . The expressions for U RR are 
easily obtained by replacing Ll A" by u R n in those for n^ 4 . As an example, we evaluate 
U RR as follows 

(3) 



U c R = E E [c'(q) i?i? M q -C'(q) i?i? M q 



q / 



with M q = Gq_ k <8> (G^A^G^), where the summation over k is restricted in the 
vicinity of the four gap nodes. For small q, M q can be evaluated to be 



127T7 2 



jq 2aro®ro-(2a-/3)ri0ri-(2a+/3)r 3 0T3 -q(Ti®T 3 + T 3 ®Ti)}, (4) 



where a = (v"j + v^)/2vfV g , (3 = (v"j — Vg)/2vfV g , and q = q • (fg + gf) with f and g the 
unity vectors parallel, respectively, to vj and v g at one of the four nodes. The upper 
and lower cutoffs of q are set to be l/l e and 1/L, respectively, where l e = ^D/2^ is 
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the elastic mean free path and L is the sample size. Substituting Eqs. (2) and (4) into 
Eq. (3), we can readily obtain IT^ = — (4/tt) ln(L// e ). Similarly, we can show that 
Ii RA = 0. As a result, we get the contribution of diagram 1(c) as 

a c s = (2/7r 2 )ln(L/Z e ). (5) 

Using the same procedure, we obtain a b a = cr|/2 = — (1/27T 2 ) ln(L/7 e ). The total cor- 
rection to the spin conductivity is given by 5a s = 2a s a + a^ + 2a s c , leading to a weak- 
antilocalization correction as 



5a s /a s = (4/a)ln(L/l e ), (6) 

where a$ = (v 2 + v 2 )/A-K 2 VfV g is the universal spin conductivity || |l8| and satisfies 
the Einstein relation <7q = pqD/4. For the electrical conductivity, each vertex in the 
diagrams equals to — evyTo|l8|], from which we can easily show that ab = a c = 0. 
Then the total correction to the electrical conductivity is given by 5a = 2a a , yielding 
a weak-localization correction as 



5a/a = -(2/a)HL/l e ), (7) 

with do = e 2 Vf/ir 2 v g the universal electrical conductivity^]. 

In the unitary limit (r/ — > 0) and at the perfect nesting (/x — > 0), there exist the 

7r-mode cooperon and diffusonfl^, p^| . Owing to the GPHS we have T2G^ A \e)T2 = 

^Q+k( e ) W1 ^h Q = (±7r/a, ±7r/a) the nesting vector. Any small deviation either from 
the unitary limit or from the perfect nesting makes the 7r-mode cooperon and diffuson 



gapped. They are given by D ff (q; e, e') RR ^ A ) = CV(q; e, e') RR ^ = C n (q; ^ ^u^^ 
Ti, where 



e, e') RR{A) = df R{A) /[Dq 2 - i(e ± e') + 25], 

with 5 = 2r/ 2 7 + fi 2 /l-f < 7 and -d RA = d RA = d RA = d RA = d RR = d RR = -d RR = 
d RR = — 47 2 /7r/?o- To the lowest order, besides the 0-mode cooperon, the 7r-mode 
cooperon also contributes to the QIE, the corresponding diagrams can be obtained by 
replacing q by Q + q in all the diagrams in Fig. 1. By summing up all the contributions 
of the 0-mode and 7r-mode cooperons, we obtain 

— = -2— = -ln 1 + — j ). 8 

At finite L, both corrections given by Eq. (8) are suppressed by decreasing 5/7 and 
vanish at 5 = 0, indicating that a s — > a$ and a — ► do. This result is in agreement 
with the numerical studies for a s in the weak-disorder limit Q]. Here we show that the 
physical origin is the existence of 7r-mode cooperon. The contributions of 0-mode and 
7r-mode cooperons have the same magnitude but opposite signs due to the GPHS. As 
a result, the corrections to the spin (electrical) conductivity coming from the 0-mode 
and 7r-mode cooperons just cancel each other out. 
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It is instructive to analyse the scattering processes described by Fig. 1. While Fig. 
1(a) yields a suppression of forward scattering of QPs, Fig. 1(b) corresponds to an 
enhancement of back scattering. Therefore, both Figs. 1(a) and 1(b) give rise to the 
weak- localization effect on the QP states. Figure 1(c) represents a more complicated 
scattering process, including an enhancement of forward scattering (k and k' located 
near the same node) and a suppression of back scattering (k and k' located near the 
oppsite nodes). The total contribution of diagram 1(c) leads to a weak-antilocalization 
effect on the QP states. As has been seen above, the QP derealization stems from the 
fact that the weak-antilocalization effect prevails over the weak-localization one. In or- 
der to understand the novel QIE, we depict in Fig. 2 the Feynman paths corresponding 
to diagrams l(a)-l(c). Figure 2(a) describes a pair of QP scattering paths, in which 




(a) (b) (c) 



Figure 2: Semiclassical scattering paths of QPs corresponding to Figs. l(a)-l(c). 

the closed loop circled twice involves only one of the two paths. Figure 2(b) represents 
a pair of paths that differ by a sequence of scattering events transversed in opposite 
directions. The scattering paths in Fig. 2(c) look like a composite of Figs. 2(a) and 
2(b). All the interference effects in Fig. 2 arise from a combination of impurity and 
Andreev scattering processes. In the unitary limit and at perfect nesting, the contri- 
butions of 0-mode and 7r-mode cooperons cancle each other. This is because the phase 
differences of coherent paths in Fig. 2 for 7r-mode cooperon differ by ir from those of 
0-mode cooperon, due to the additional GPHS. 

Since the QP spin is conserved, the behavior of the dimensionless spin conduc- 
tance g s = a s /(l/2) 2 as a function of L enables construction of a scaling theory of 
QP (derealization]!]. The one-parameter scaling function is defined by (3{g s ) = 
din g s /din L. It is easily shown that both Eqs. (6) and (8) can be collapsed into 
a single universal scaling function as 

P(g s ) = 8/vrV , (9) 

for g s — > oo. The positive (3{g s ) strongly indicates the existence of extended low-energy 
QP states. This spin metal state is characterized by the absence of charge diffusion, 
as the electrical conductivity has a weak-localization correction. The observation that 
(3{g s ) decreases with g s implys that these extended low-energy QP states are different 
in character with the usual extended bands. We argue that this novel phenomenon 
is related to the strong anisotropy of energy gap in cf-wave superconductors, as the 
wave vectors of all extended low-energy QP states are nearly along the four nodal 
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directions. It is worthy to mention that the extended low-energy QP states have been 
predicted by Balatsky et from the novel network of delocalized impurity states 

in the unitary limit. However, the physical mechanisms for the formation of extended 
states are different from each other in Ref.(7| and in the present theory. On the other 



hand, the nonlinear-sigma-model calculations in Refs.|l], [|, 12] are incomplete in which 
the contribution of Fig. 1(c) was omitted (see Fig. 1 in Ref.jj^]), and thus yielded 
a negative logarithmic correction to o~ s for symmetry class CI. The novel quantum 
interference process is expected to exist in superconductors that belong to symmetry 
classes C and D. We note that RR-cooprons are not influenced by the time-reversal 
breaking Q |3], |l2|. How a magnetic field (or dilute magnetic impurities) affects the QP 
derealization effect is another interesting and open problem. 

The QP derealization effect is expected to have a manifestation in the low-temperature 
thermal transport property. Since the QP energy is also conserved, the electronic ther- 
mal conductivity K should obey the Wiedemann- Franz law K/Ta s = Ait 2 k 2 B /?>. Here 
the temperature dependence of the quantum interference correction 5cr s (T) results from 
the dephasing time 70 (T), the latter may be obtained by considering the interactions 
between QPs[13]. Since the self-consistent T-matrix approximation is valid only for the 
case of dilute impurities^], we do not rule out the possibility of localized low-energy 
QP states at higher impurity concentrations. Should it appear, there might exist a 
quantum transition from spin metal to spin insulator in 2D disordered d-wave super- 
conductors, which is expected to be observed in a low-temperature thermal transport 
experiment. By increasing the impurity concentration, the temperature dependence of 
K would change from metallic to insulating behavior. 

In conclusion, we have studied the QIEs on the low-energy QP transport in weakly- 
disordered 2D cf-wave superconductors both in Born and in unitary limits. We find 
a new, one- loop diagram related to the LPHS, which qualitatively modifies the usual 
weak-localization results. In generic situations, the existence of this novel quantum 
interference is shown to yield a weak-antilocalization correction to the spin conductiv- 
ity. With coming close to unitarity and the nesting, this correction is supressed and 
eventually vanishs due to the GPHS. With a universal one-parameter scaling function, 
we show that the novel QIE can produce the extended low-energy QP states in the 
weak-disorder limit. 
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